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1 Introduction 



Rare decays of beauty mesons (B ± and B°) are sensitive probes of new particles which arise in 
models beyond the Standard Model (SM). In the context of supersymmetric (SUSY) extensions 
of the SM, the measurements of decays such as b — > 57, B ± — > tv, and B ± — > Dtv provide 
important constraints on the masses of new particles which are too heavy to be produced 
directly. Some recent analyses showing constraints on the parameter space of the Minimal 
Supersymmetric Standard Model (MSSM) can be found in [Ifj6]. Of much interest for the 
LHC experiments is the unobserved decay B s — > ji + n~ . Due its distinct signature, this decay 
can be searched for by three LHC collaborations: LHCb, CMS and ATLAS. As pointed out 
in [7|j9], B s — > is a very effective probe of SUSY models with large (> 30) tan/3, and its 

importance has been emphasised in numerous studies. The upper limit on the branching ratio 
(BR) of B s — > has been steadily reduced during Run II at the Fermilab Tevatron. As of 

the year 2010, limits of the order of BK(B S — > /U + /i~) < 5 x 10~ 8 (i.e. an order of magnitude 
above the prediction of the SM) were obtained by both the CDF [10] and DO 11 collaborations. 



Recently, the CDF collaboration announced a possible first signal 12 , although with a 
low significance. This result has not been confirmed by the recent searches at LHCb [13) and 
CMS (l4j. These improved limits for BK(B S — > /U + /i~) further constrain the SUSY parameter 
space, and we show in this paper that such constraints can be superior to those which are 
obtained from direct searches for squarks and gluinos. Using a combination of the individual 



limits on BR(_B S — > ) from LHCb and CMS 15 , we present an updated study of the 



constraints in the context of five distinct SUSY models. Our numerical analysis is performed 



with Superlso v3.2 16-18 , and we study the following SUSY models: the constrained MSSM 
(CMSSM), non-universal Higgs mass (NUHM), anomaly mediated supersymmetry breaking 
(AMSB) and gauge mediated supersymmetry breaking (GMSB), all in the context of the MSSM; 
we also study a semi-constrained version of next-to-MSSM (NMSSM). Moreover, we consider 
an alternative observable which includes BK(B S — > /i + /i~), namely, a double ratio of leptonic 



decays 19,20 . The double ratio has no dependence on the absolute value of the decay constant 



Jb s , which is the main source of uncertainty in the SM prediction for BR(_B S — > and it 

was shown in 21 that this observable can provide competitive (or even superior) constraints 
on the SUSY parameter space. The main uncertainty in the SM prediction for the double ratio 
arises from the CKM matrix element \V u b\, for which the prospects of precise measurements at 
high-luminosity B factories are very promising. The final integrated luminosity of the operation 
of the LHC at = 7 TeV is likely to be significantly larger than the amount which was 
anticipated at the start of the run, which could enable the SM prediction for BK(B S — > /i + /i~) 
to be probed. We discuss the expected sensitivity to BR(B S — > /U + /i - ) as a function of the 
integrated luminosity, as well as the prospects for a measurement of a SM-like BK(B S — > /i + /i~) 
during the 7 TeV run, and its potential implications. 

Our work is organised as follows: in sections 2 and 3 we present a theoretical introduction 
to the decay B s — )■ and the double ratio respectively; section 4 contains our numerical 

analysis of the constraints on various SUSY models that are obtained from the recent upper 
limit on BR(B S — > /i + /i~); in section 5 the experimental prospects for B s — > ji~ are discussed, 
and conclusions are contained in section 6. 
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2 The decay B s — > /i + \i 



It has been emphasised in many works [l]H[7l|^[22}|24] that the decay B s —> fi + fi is very 
sensitive to the presence of SUSY particles. At large tan/3, the SUSY contribution to this 
process is dominated by the exchange of neutral Higgs bosons, and very restrictive constraints 
are obtained on the supersymmetric parameters. The BR(B S — > can be expressed 

as fT7l[25H27 



BR(B S /iV 



647T 3 
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Ami 



m 
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where the coefficients Cq 1 , Cq 2 , and Cio parametrize different contributions. Within the SM, 
Cq 1 and Cq 2 are negligibly small, whereas the main contribution entering through C\o is helicity 
suppressed. In SUSY, both Cq 1 and Cq 2 can receive large contributions from scalar exchange, 
which was first pointed out (in the context of a different decay, b — > sl + l~) in 28 . The explicit 
expressions for the different coefficients can be found in e.g. [TF| . 

The B s decay constant, fB s , constitutes the main source of uncertainty in BR(_B S — > 
As of the year 2009 there were two unquenched lattice QCD calculations of fs a) by the HPQCD 
and FNAL/MILC 

Ml 



29 



collaboration 
f Bs = 238.8 ± 9.5 MeV 



[30] respectively, which when averaged gave the value 
The calculation of |30| was updated in 32 , which gave rise to a 
higher world average of Jb s = 250 ± 12 MeV in the year 2010. Recently, the ETM collaboration 
announced its result of fs a = 232 ± 10 MeV [33]. At the Lattice 2011 conference [34], new 
results by FNAL/MILC (f Bs = 242 ± 9 MeV) and the HPQCD collaboration (f Bs = 226 ± 10 
MeV [35] and f Bs = 225 ± 4 MeV |36|) suggest that an updated world average would be lower 

238.8 ± 9.5 MeV 



31 



than that of the year 2009. In our numerical analysis we will use f Bg 

To study the constraints on the parameter spaces of SUSY scenarios, we use the newly 
released combined limit from LHCb and CMS at 95% C.L. [15]: 



BR(B S /U + /0 < 1.1 x 10" 



(2) 



More details are given in section [5] In order to take into account the theoretical uncertainties, 
in our numerical analysis we will use the following limit 



BR(B S < 1.26 x 10~ 8 

to constrain the parameter spaces of the SUSY models under consideration. 



(3) 



3 The double ratios of purely leptonic decays 

The main uncertainty in the theoretical prediction of B s — > fi + n~ is f Bs . As described in 
Section |2j f Bs is now being evaluated in the unquenched approximation by various lattice 
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collaborations. The error (which is currently around 5% or less) has been reduced over time, 
and the central values of fs a from the various collaborations are in reasonable agreement. 
The prospects for a further reduction of the error are good. However, despite the continuing 
improvement in the calculations of fs a our view is that it is instructive to consider other 
observables which involve B s — > ji + n~ but do not depend on the decay constants, and to 
compare the constraints on the SUSY parameter space with those which are obtained from 
BR(£> S — > alone. One such observable which involves B s — > but has essentially 

no dependence on the absolute values of the decay constants, is a double ratio involving the 



leptonic decays B u — > tu, B s — > , D — > fiu and D s — > hu/tv 19,20 
One such double ratio is defined by: 

T(B S -> n+p-) T(D -)■ H \V ts V tb \ 2 a 2 (f D /f Ds ] 



Y(B U ^TU) T(D s ^fIls) \V ub \ 2 7T 2 (f B /fB a ) 2 



(4) 



The quantity (fs / /b s )/(/d/ /d s ) deviates from unity by small corrections of the form m s jm b 
and m s /m c . The double ratio would be equal to one in the heavy quark limit of a very large 
mass for the b and c quarks, and in the limit of exact SU(3) flavour symmetry (m s — > 0). 



A calculation in 19 gives (/b//b s )/(/d//d s ) = 0.967, and subsequent works [37] also give 
values very close to 1, with a very small error. Unquenched lattice calculations of the ratios 
/d s //d and JbJ }b have a precision of the order of 1% (e.g. [32]), from which it can be inferred 
that the numerical value of the double ratio is very close to 1. In our numerical analysis we 
will take (/#/ '/b„)/(/d/ Jd s ) = 1- Importantly, the absolute values of the decay constants 
do not determine the value of the double ratio, in contrast to the case of BR(£> S — > /i + /i~) 
alone in Eq. ([I]). Instead, \V u b\ replaces fs s as the only major source of uncertainty, as can be 
seen from Eq. Q. Information on \V u b\ is available from direct measurements of semileptonic 
decays of B mesons, both inclusive (B — > X u tv) and exclusive (B — ^ -kIv). Moroever, global 



fits 38 in the context of the SM give additional experimental information on \y u \\- Due to its 
different theoretical uncertainties, the double ratio is an alternative observable which includes 
B s — > and can provide competitive constraints on SUSY parameters. A comparison 

of the constraints on specific SUSY models from the double ratio and from BR(i? s — > /i + /i _ ) 
alone is of much interest because the theoretical input parameters \V u b\ and fs s for these 
two observables are independent. Such a comparative study was performed for the first time 
in [2l], and it was shown that the double ratio can provide stronger constraints. In particular, 
the constraints from the double ratio are maximised (minimised) for smaller (larger) \V u b\, while 
the constraints from BR,(B S — > /i + /U~) are maximised (minimised) for larger (smaller) fs s - 
We will perform an updated study of these two observables using the recently improved 



bounds on BR(5 S — > ) from the LHCb [13] and CMS 14 collaborations. In the previous 



study of the double ratio plj the value \V ub \ = (3.92 ± 0.45 ± 0.09) x 10" 3 39 was used 



which is an average of the exclusive and inclusive determinations of \V u b\. We note that this 
world average does not include three recent measurements of \V u b\, of which two are from the 



exclusive channel 40 41 and one is from the inclusive channel 42 . The inclusion of these 
measurements would only have a small effect on the world average, and so for simplicity we 
will use I Kb I = (3-92 ± 0.45 ± 0.09) x 10~ 3 , as done in [2l]. 



We note that the exclusive determination of \V u b\ suggests values of \V u b\ which are below 
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the central value of the world average. The exclusive determination of \V u b\ requires a theo- 
retical calculation of one hadronic form factor f+(q 2 ) (where q is the momentum of £). For 
q 2 > 16 GeV 2 one can use lattice QCD to calculate f + (q 2 ), while for q 2 < 16 GeV 2 non-lattice 
techniques must be used. In both regions of q 2 the extracted value of \V u b\ is below the central 
value of the world average. The inclusive determination of \V u b\, which does not have a depen- 
dence on lattice QCD, suggests values which are above the central value of the world average. 
Prospects for precise measurements of \V u b\ in the inclusive channel are very good at high- 
luminosity B factories. In particular, the method used in 42 43 is a very promising approach 
because the theoretical errors are greatly reduced by employing a low cut on the momentum 
of the t (j)£ > 1 GeV), which keeps 90% of the phase space of B — > X u £v. This anticipated 
experimental improvement in the measurement of \V u b\ bodes well for the double ratio as an 
alternative observable with which to constrain SUSY. It is important to emphasise that fs s is 
currently known with greater precision than \V U b\, and this may also be the case in the era of 
a high-luminosity B factory. However, we note that the central values of these unrelated input 
parameters plays a major role in determining which observable gives the stronger constraints, 
as will be discussed in our numerical analysis. 

The double ratio also has the attractive feature of using ongoing measurements of BR(D S — >■ 
Hv/tv) and BR(D — > /iz/). Such decays are not usually discussed when constraining SUSY 
parameters (although see [44] for a discussion of D s — > jiv/ru in this regard), but increased 
precision in their measurements would enhance the capability of the double ratio to probe 
the SUSY parameter space. The decay B u — > tv alone is very sensitive to the presence of a 
charged Higgs boson (H ) and provides a strong constraint on tan (5 and the mass of H ± in 
SUSY models 45 -48]. The experimental prospects for precise measurements of all the decays in 
the double ratio are very promising. The precision in the measurements of BR(D S — > fiv) and 
BR(_D S — > tv) will be improved at the ongoing BES-III experiment [49], and at high-luminosity 
B factories operating at a centre-of-mass energy of y/s ~ 10.6 GeV (and also possibly at energies 
in the charm threshold region). Similar comments apply to the prospects for significantly 
improved measurements of BR(D — > fiv) and BK(B U — > tv). For more details about the 
calculation of these decays we refer the reader to [21] . 



In this analysis we use 

As- 5 -, (5) 

VSM 

where 

/ BR(B fl ->//+//- ) \ / / BR(£>g tv) \ 
71 ~ V BR(B U ^tv) ) I V BR(£> ^[iv))' U 

The theoretical evaluation of //sm gi ves (2.47 ± 0.58) x 10~ 7 where the main uncertainty 
comes from V u b- To determine the experimental limit on the ratio R, we combine the limits 
on the individual branching fractions, namely BK(B S — > yU + /i _ ), BK(B U — > tv), BK(D s — > tv) 
and BR(D — > jj,v). To compute the p. d.f of R, we use a Gaussian distribution for the measured 
decays, and a "truncated" Gaussian p. d.f for the upper limit in We consider two different 
approaches. The first approach consists in building first the p.d.fior BK(B S — > fi + ^r) which 
reproduces the 90% and 95% C.L. experimental limits, and to combine it with the Gaussian 
p. d.f of the other involved decays. The second approach determines the p. d.f 'from the derivative 
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in 0.12 




Figure 1: Probability distribution function (p.d.f) of the double ratio R. For B s — > the 
p. d.f was obtained based on the CL. from [15]. For the other three decays, their measurements 
are modelled as Gaussians. 



of the CL s+ b with BR(S S — > n + fi ), which is extracted from the derivative of CL S shown in 
Ref. 15 and the almost constant behaviour of CL^ . Fig. [T] shows the R p.d.f . 



Both approaches agree and provide the upper limit for R, at 95% C.L.: 

R < 2.3 , (7) 

in which the uncertainty from V u b is taken into account. In our numerical analysis we use ^ to 
constrain the supersymmetric parameter space in various scenarios in the MSSM and NMSSM. 



4 Constraints on SUSY Models 



We consider five distinct SUSY models in order to illustrate the impact of the new limits on 
BR(B S — > on the SUSY parameter spaces. All previous studies have been carried out 

limits were released^ Some very recent works 51 study 



before the LHCb 1 13 1 and CMS 



14 



the impact of the latest CDF result 12 only, and address the case of the excess of events being 



a genuine signal. Moreover, none of the previous studies have considered the double ratio, apart 



from our earlier work in 21 in which two of the five SUSY scenarios were discussed. 



For each scenario we also check the constraints from direct searches for Higgs bosons and 
delimit the regions where the lightest supersymmetric particle (LSP) is charged. All the flavour 



observables are calculated with the Superlso v3.2 program 16418 . The spectrum of the MSSM 



points is generated with SOFTSUSY-3.1.7 |52| and we used NMSPEC program from the 
NMSSMTools 3.0.0 package [53] for the NMSSM points. For every generated MSSM point 



we check if it fulfills the constraints from the Higgs searches using HiggsBounds-3.2.0 54 55 
The value of m t = 173.3 GeV [56] is used throughout. 



In an updated version of Ref. 50 the impact of the latest LHCb FLs] and CMS [l4] limits on BK(B S — s- ) 



is studied amongst other observables in a global CMSSM fit. 
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4.1 CMSSM 



The first model we consider is the constrained MSSM (CMSSM) [57], which is characterized 
by the set of parameters {m , my 2, A , tan/3, sgn(/x)}. The CMSSM model invokes unifica- 
tion boundary conditions at a very high scale vtigut where the universal mass parameters are 
specified. 

To explore the CMSSM parameter space, we generate about 300,000 random points scanning 
over the ranges m G [50,2000] GeV, m 1/2 G [50,2000] GeV, A G [-2000,2000] GeV and 
tan/3 G [1, 60] with positive \i (as favoured by the muon (g — 2) measurements). 

The results are displayed in Fig. |2j where the four-dimensional space is projected into a 
plane. When interpreting these results it is therefore important to remember that each point 
in the figures corresponds to a multi-dimensional parameter space in the variables which are 
not displayed on the x-axis and the y-axis. 

In order to show the viable parameter space of the SUSY scenario under investigation, 
in all the figures we introduce a colour coding which is applied sequentially. Areas which are 
disallowed theoretically are in white. Next, the points which are disallowed phenomenologically 
are plotted, which are those with a charged LSP (in violet) and those which are excluded by the 
direct searches for Higgs bosons (in black). In this way, these points lie in the background. On 
top of them, the points excluded by the double ratio R (in orange) are displayed, superseded 
by the points excluded by BR(5 S — » /x + ji^) (in yellow). Finally the allowed points (in green) 
are shown in the foreground. 

These indirect constraints on the CMSSM parameter space from BR(5 S — > yU + /i~) are 
competitive with the direct constraints from searches for squarks and gluinos by ATLAS and 
CMS [58] . As expected, one can see strong constraints on small and large tan j3 values. At 
large tan/3 (> 30), these constraints are stronger than those obtained from BR(£> — > X s j) [3]. 

In order to better quantify the impact of BR(B S — > /i + /i~) and R, we show in Fig. [3] the 
constraints for fixed values of tan /3 (=30, 40 and 50) and Aq = 0. One striking result here 
is that the double ratio, being a combination of four different flavour observables, extends 



impressively the constraints obtained by BR(i? s — > ji + jj, ) alone, as was pointed out in 21 



Also, for tan (3 = 50, the constraints from the flavour observables go far beyond the direct 
search limits by the ATLAS and CMS collaborations for the same scenario. 

The SUSY contributions to B u — > tv gives rise to a scale factor which multiplies BK(B U — >■ 
tu). When we manually set this scale factor to be equal to 1 (as in the SM), the excluded 
region of the plane [mo, m^] does not change much. Therefore we conclude that the points i) 
and ii) above are the main reasons why the double ratio gives the superior constraints. 



4.2 NUHM 

The second model we consider involves non-universal Higgs masses (NUHM) [59]. This model 
generalizes the CMSSM, allowing for the GUT scale mass parameters of the Higgs doublets 
to have values different from m , i.e. ^ ^h 2 m o- These two additional parameters 

with dimension of mass can be traded for two other parameters at a lower scale, which can be 
conveniently chosen as the \x parameter and the mass rriA of the CP-odd Higgs boson. 
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Figure 2: Constraints from BR(i? s — > and the double ratio R in the CMSSM planes 

(mi/2, mo) in the upper panel, (rrig, rriu L ) in the middle panel and (m^, tan (3) in the lower panel. 
The colour coding is given in the text and the constraints are applied in the order they appear 
in the legend, with the allowed points in green displayed on top. 
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Figure 3: Constraints from BR(5 S — > fi + fi~) and the double ratio R in the CMSSM planes 
(mi/ 2 ,m ) on the left and (m^m,^) on the right, for A = and tan/3 = 50 (upper panel), 
tan/3 = 40 (middle panel) and tan/3 = 30 (lower panel). 
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Figure 4: Constraints from BK(B S — > ) and the double ratio R in the NUHM plane 



(rattan /3). The colour coding is explained in section 4.1 




Figure 5: Constraints from BK(B S — > and the double ratio R in the NUHM parameter 

plane (/i, m^) with tan/3 = 35, A = 0, m = 1000 and mi/ 2 = 500 GeV on the left, and in the 
plane (m#+,tan/3) for fi = 1000, A = 0, m = 800 and = 600 GeV on the right. 



We generate about 300,000 random points in the ranges m G [50,2000] GeV, G 
[50,2000] GeV, A e [-2000,2000] GeV, tan/3 G [1,60], /i G [-2000,2000] GeV and m A G 
[20, 1000] GeV. The results are presented in Fig. |4j Again the constraints are very important, 
and restrict strongly the region of large tan/3 / small m^. 

In Fig. [5] we show two examples in the two-dimensional parameter planes m^) and 
(m,H+ , tan /3) with the rest of parameters being fixed. As can be seen from the figures, a large 
part of the parameter space is restricted by BK(B S — > yU + /i~) and R observables, whereas in 
the same plane one would not get any constraints from BR(i? — > V s 7) for fj, > (3|. 
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Figure 6: Constraints from BR(5 S — > fi + fi ) and the double ratio R in AMSB. The colour 



coding is given in section 4.1 




Figure 7: Constraints from BR(5 S — > fi + fi ) and the double ratio R in the AMSB parameter 
plane (m 3/ / 2 ,tan/3) for m = 500 GeV and ji > (on the left) and \i < (on the right). 



4.3 AMSB 

We can now focus on another supersymmetry breaking scenario, namely the Anomaly Mediated 
Supersymmetry Breaking (AMSB) [60j. This is a special case of gravity mediation in which 
there is no direct tree-level coupling that transmits the SUSY breaking in the hidden sector 
to the visible one. The breaking is communicated through the conformal anomaly. The free 
parameters of the minimal model consist of {m , m 3 / 2 , tan/3, sgn(/i)}. 

Previous studies were performed in [4|[6l] . To explore the parameter space of AMSB, we 
scan over m G [50,2000] GeV, m 3/2 e [1,100] TeV and tan/3 e [1,60], and generate 300,000 
random model points. The results are presented in Fig. [6] and show stronger limits for low 
values of mo and large tan (3. 
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Figure 8: Constraints from BR(/3 S — > fi + fi ) and the double ratio R in the AMSB parameter 
plane (m ,tan/3) for m 3 / 2 = 30 TeV and /i > (on the left) and /i < (on the right). 



Fig. [7] and Fig. [8] show the results in two-dimensional planes in order to see better the extent 
of the constraints. In the plane (7773/2, tan /3) for mo = 500 GeV, essentially all of the points 
with tan/3 > 30 are disfavoured regardless of the value of 7773/2. In the plane (m ,tan/3) with 
7773/2 = 30 TeV, one obtains strong constraints for small / large tan /3. Scenarios with \i < 
show similar effects, with the constraints being less pronounced. It is also evident that a large 
portion of the parameter space is already excluded by the constraints from the direct searches 
for Higgs bosons, as implemented in the HiggsBounds program. 



4.4 GMSB 

The last MSSM scenario that we consider is the Gauge Mediated Supersymmetry Breaking 
(GMSB) scenario [62), which consists of the SUSY breaking sector and the messenger sector. 
The latter can be taken as a 5 + 5 of the SU(5) which contains the Standard Model group, and 
therefore the gauge coupling unification is not affected. The minimal model is characterized by 
the set of parameters {A, M mess , N 5 , c grav , tan/3, sgn(/i)}. For our study, we consider N 5 = 1, 
Cgrav = 1 and generate about 300,000 random points in the ranges A e [10, 500] TeV, M mess G 
[10 2 , 10 14 ] TeV and tan /3 e [1, 60] with A < M mess . 

In Fig. [9] we show the results in the parameter planes (A, tan/3) and (M mess , tan /3). Again, 
the region of large tan /3 is the most restricted by the flavour observables. To see better the 
regions in the parameter space which are excluded by BK(B S — > /ii + /U~) and the double ratio 



.R, we show in Fig. 10 the results in the plane (M mess , tan /3) for a fixed value of A = 100 TeV 
for both fM > and /i < 0. It is remarkable to see that tan (3 > 40 is excluded regardless of the 
value of M mess , while the same plane is probed by the well-known BR(73 — > X s ^) constraints 
only for a very large messenger scale (M mess >10 10 TeV) |3|. 
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Figure 9: Constraints from BR(i? s — > and the double ratio R in the GMSB parameter 

planes (A,tan/3) on the left and (M mess , tan (3) on the right. The colour coding is given in 
section I4.ll 




Figure 10: Constraints from BR(S S — > fi + fi ) and the double ratio R in the GMSB parameter 
plane (M mess , tan /3) with A = 100 TeV, for fi > (on the left) and [i < (on the right). 



Fig. 11 presents the constraints in the plane (A,tan/3) with M mess = 500 TeV and shows 
that only relatively small values of A are affected by B s — > fi + fi~. In the white area (which is 
especially large in the case of n < 0) it is not possible to find any valid model point. 
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Figure 11: Constraints from BR(5 S — > fi + fi ) and the double ratio R in the GMSB parameter 
plane (A,tan/3) with M mess = 500 TeV, for /i > (on the left) and \i < (on the right). 



4.5 CNMSSM 



The last scenario that we consider is a constrained version of the NMSSM (CNMSSM) with 
semi-universal parameters defined at the GUT scale 63 . The choice of a semi-universal scenario 



instead of the case of strict universality facilitates the obtention of valid NMSSM points [64 j. 
In this scenario, k, A and m| are computed from the minimization equations and the free 



parameters are {m , m^, A 0} A K , A, tan j3, sgn (//)}. Previous studies were performed in 65 66 
Our sample of 300,000 random points is generated in the ranges m G [50, 2000] GeV, mi/2 G 
[50,2000] GeV, A G [-2000,2000] GeV, A K G [-2000,2000] GeV, A G [10~ 3 , 1] and tan/3 G 
[1,60]. 



The results are displayed in Fig. 12 in the parameter planes (m#+,tan/3) and (A,tan/3). 
The constraints are more severe for large tan /3, small m#+ and large A. In Fig. 13 we fix two 
of the parameters, namely A = 0.01 and tan /3 = 50. This allows us to see in a clearer way the 



effect of the constraints on the other parameters. In Fig. 14 the same results are shown for 
A = 0.1. As mentioned before, the constraints are more pronounced for larger A. 

As a final example we fix all the parameters except two, to see the results in a two- 
dimensional plane. This is done in Fig. 15 for A = 1000 GeV, A K = —60 GeV, tan j3 = 50 and 
A = 0.1. As can be seen, a large part of this parameter plane is excluded by BR(_B S — » fi + fi~) 
and the double ratio R. 
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Figure 12: Constraints from BR(i? s — > fi + fi ) and the double ratio R in CNMSSM in the 



parameter planes (mff+, tan f3) and (A,tan/3). The colour coding is given in section 4.1 



4.6 Discussion 

In the above subsections, we investigated the constraining power of BR(B S — > fi + fi~) and the 
double ratio R for different SUSY scenarios. As explained in sections [2] and |3j the main input 
parameter for BR(B S — > H + is fs 3 , while \V u b\ is the most important input for R. To 
examine how the choice of these input parameters can affect our results, we consider here two 
scenarios, namely the "least constraining" (with high \ V u b\ and low fs a ) and "most constraining" 
(with low 1 1^6 1 an d high f Bs ) cases for BR(£> S — > and the double ratio. For the least 

constraining scenario we consider the inclusive determination of \V u b\ with the central value 
being 4.34 x 10~ 3 39 , and fs a = 232 MeV 33 . For the most constraining case we take 



the exclusive value |T4&| = 3.42 x 10~ 3 39 and fg a = 250 MeV 32 . To compare these two 
cases we take an example in the CMSSM scenario with tan/3 = 40 and A = 0. The results 
are presented in Fig. [16J As can be seen, in the most constraining case, the exclusion limits 
are greatly increased while in the least constraining case the results are only slightly changed. 
This shows that the analysis in the previous subsections does not correspond to a particularly 
optimistic choice of the input parameters. 

The next point we discuss here is the effect of B u — > rv in the double ratio. The constraints 
from B u — > rv alone on the parameter space of [m , my 2 ] have been presented in several works 
(e.g. [17]) and the excluded region differs from that obtained from BR(_B S — > alone, 



as can be seen in Fig. 17 In most of the parameter space, BR(B U — > rv) is reduced with 
respect to the SM value, leading to the large blue excluded strip in Fig. [TTj On the other hand, 
in the small strip, BR(i?„ — > rv) is larger than in the SM. In the narrow region in between, 
a cancellation happens since the charged Higgs contribution is roughly twice that of the SM 
contribution and so B u — > rv cannot exclude this parameter space. As can be seen from the 
figure, BR(£> S — > n~) probes larger values of mi/ 2 than B u — > rv, although B u — > rv can 
exclude part of the region 1300 GeV < m < 1600 GeV and m^j < 200 GeV which cannot be 
excluded from BR(_B S — > //+ alone and the double ratio. 
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Figure 13: Constraints from BK(B S — > fi + {i ) 
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Figure 14: Constraints from BK(B S — > fi + fi ) and the double ratio R in different CNMSSM 
parameter planes for A = 0.1 and tan/3 = 50 with jj, > 0. 
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Figure 15: Constraints from BR(5 S — > and the double ratio R in the CNMSSM param- 

eter plane (m 1 / 2 ,fn ) for A = 1000 GeV, A K = —60 GeV, tan/3 = 50 and A = 0.1. 
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Figure 16: Constraints from BR(5 S — > and the double ratio R in the CMSSM parameter 

planes (my 2 ,mo) f° r tan/3 = 40. On the left, the most constraining case with low \V u f,\ and 
high Jb s and on the right the least constraining case with high \V u b\ and low /# s . 



The reason why the double ratio is more constraining than BK(B S — > alone is mainly 

due to two reasons: i) \V u b\ is used as an input parameter in the double ratio, instead of /s s . 
Although these two parameters have comparable errors, their current central values give rise 
to stronger constraints from the double ratio, as discussed in the preceding paragraph. This 
could not have been expected, and a value of fs s much larger than that preferred by lattice 
QCD would have ensured that BR(5 S — > alone had the stronger constraints; ii) The 

experimental value of BR(I? U — > rv\ which enters the derivation of i] in Eq. ([6]), is larger 
than the SM expectation, and so reduces R in Eq. ([7]), leading to a stronger constraint on the 
SUSY parameter space. The SUSY contributions to B u — > tv gives rise to a scale factor which 
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Figure 17: Constraints from BR(5 S — > the double ratio R, and B u — > tv in the CMSSM 

plane (m,i/ 2 , mo), for tan/3 = 50 and A = GeV. This figure is the same as Fig. [2] (upper left 
panel) but with the constraint from B u — > tv superimposed. 



multiplies BR,(B U — > tv). When we manually set this scale factor to be equal to 1 (as in the 
SM), the excluded region of the plane [m , my^ does not change much. Therefore we conclude 
that the points i) and ii) above are the main reasons why the double ratio gives the superior 
constraints. 

Finally we discuss the effect of a hypothetical measurement of BR(5 S — > fi + ^) at the SM 
value (3.5 ± 0.3) x 1CT 9 . Fig. 18 shows the obtained impact in the CMSSM plane (m^tan/?) 
with all the parameters being varied in the intervals given in section 4.1 For comparison, the 



same parameter plane with the current experimental limits is also provided. As can be seen, 
almost no scenario with tan/3 > 45 remains viable regardless of the other parameters in the 
case of a SM-like discovery, and the parameter space of the CMSSM becomes very restricted. 



5 Experimental prospects 

At present, the best upper limit for BR(5 S — > /i + /U~) measured in a single experiment comes 

(8) 



fromLHCb 13 



BR(B S /i + /i~) < 1-5 x 10" 



at 95% C.L. This upper limit is followed closely by the result from CMS 14 

BR(B S ->• < 1.9 x 10" 8 



(9) 



at 95% C.L. These two results were officially combined for EPS conference in Ref. 15 , giving 
the upper limit of 

BR(B S < 1.1 x HT 8 , (10) 
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Figure 18: Constraints from BR(B S — > /i + /i~) and the double ratio i? in the CMSSM parameter 
planes (m^ , tan /3) in the hypothetical case of a SM-like measurement (lower panel) and with 
the current experimental limits (upper panel). In the left panel the allowed points in green are 
displayed in the background while in the right panel they are in foreground. 



which we will use to constrain the parameter space of SUSY models. The CDF collaboration 
obtains a 95% C.L. upper limit [i~2l: 



BR(B S -»■ [i+n-) < 4.0 x 10~ 8 , (11) 
together with a one sigma interval 

BR(B S -> fi+fT) = (1.8+li) x 10" 8 , (12) 
coming from an observed excess over the expected background which corresponds to a p— value 



of 0.27%. Finally, the DO collaboration obtains the 95% C.L. upper limit 11 : 

BR(B S -> < 5.1 x KT 8 . (13) 

The preliminary result on BR,(B S — > /i + yU~) |15| from the combination of the limits from 
LHCb and CMS shows an excess of more than one sigma (CL b w 0.92 for values of the BR 
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around the SM value) with respect to the background-only hypothesis. This excess can be 
accounted for by a BR(5 S — > ~ (3.7^' 7) x 10 -9 . However, the signal significance is not 

enough to claim evidence. In this section we study the experimental sensitivity to B s — > [i + 
and the prospects for its measurement in the period of operation of the LHC at ^/s = 7 TeV. 

5.1 Combination LHC-CDF 

The CDF experiment at the Tevatron has reported a p— value of 0.27% for the background only 



hypothesis 12 . In order to evaluate whether a combination of results from CMS, LHCb and 
CDF could lead to evidence for a signal, we perform an approximate combination of the results 
of the three experiments, based on the signal and background expectations and the observed 



pattern of events. We use mc_limit 67 to combine the results of the different experiments and 



to extract the confidence levels. We have also scaled fd/f s to the value measured at LHCb 68 
in order to be consistent with the value used in the LHC combined result. 

According to this study, a hypothetical combination of the LHCb and CMS results with that 
of CDF would increase CL^ to ~ 0.994 (for values of the BR close to the most probable value), 
which is close to a 3a deviation. Note that this is approximately the same signal significance 
that CDF obtains alone. This approximate study leads to the following averaged branching 
ratio: 

BR(B S ->• i2 + fi-) C DF + LHC w (6±|) x 10- 9 . (14) 

However, at the time of writing this paper, this kind of combination has not been performed 
officially 

5.2 Sensitivity to BR(B S -> /j+fi' ) at the LHC 

We perform a toy MC study in order to determine how much luminosity is needed to obtain 
evidence for B s — > at the LHC. For this, we scale the signal and background expectations 

accordingly with the increase of luminosity. Fig. [19] shows the integrated luminosity that is 
needed in order to obtain a 3(5) a evidence (discovery) of a given BK(B S — > /U + /i~) in either 
LHCb or CMS. 

Assuming that the ratio of luminosities between CMS and LHCb remains at the value of 
the current analysis (i.e. CMS takes approximately four times more data than LHCb over 



the same period of time), we show in Fig. 20 the integrated luminosity scale factor (with 



respect to the amount of data used in [15]) that would be needed for the discovery of a given 
BK(B S — > /i + /i~) in the case of a CMS+LHCb combination. The width of the bands reflects 
possible scenarios for the evolution of the systematic uncertainties, where the lower side assumes 
negligible systematics and the upper side assumes that the systematics do not get reduced 
with time. It can be seen that with 6-8 times more luminosity than that used in Ref. [15] 
a CMS+LHCb combination could provide evidence at the 3a level for BK(B S — > yU + yU~) of 
the SM. This corresponds to between 2 and 3 fb _1 for LHCb and between 7 and 10 fb _1 for 
CMS. As the sensitivity of CMS is equivalent to that of LHCb for four times more luminosity, 
a scenario in which CMS takes up to 14 fb _1 and LHCb takes 2 fb -1 would afford equal 
sensitivity as a combination of CMS with 10 fb" 1 and LHCb with 3 fb _1 . From this toy MC 
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Figure 19: Required luminosity in order to provide a 3a evidence (orange) or a 5a discovery 
(green) of a given BR,(B S — > on the left for LHCb and on the right for CMS. 



study we conclude that the SM prediction for BR(B S — > is likely to be probed during 

the operation of the LHC at = 7 TeV (i.e. before the end of the year 2012). If ATLAS 
can manage to obtain sensitivity to BR(5 S — > /U + /i~) which is comparable to that of CMS, 
then even a 5a discovery for a SM-like BR(I? S — > fi + n~) would be possible during the run at 



a/s = 7 TeV. However, from pre-LHC MC studies in Ref. 69 the sensitivity of ATLAS was 
found to be inferior to that of CMS. If experimental evidence of B s — > is achieved at the 

LHC, the double ratio in Eq. Q would be measured for the first time. Moreover, limits on the 
ratio BR(i?d — > fi + n~)/BR(B s — > /U + /i~) (which is a very interesting test of Minimal Flavour 
Violation) would also be set. If BK(B S — > / u + ^r) is much smaller than the SM prediction (as 
can happen for example in the MSSM [70] and NMSSM), values down to 0(5 x 1CT 10 ) can still 
be discovered with an upgrade of the LHCb. 

5.3 NP discovery with B s — > 



In section |5.2 we discussed the luminosity needed for discovery of B s — > fi + . However, a 
measurement of B s — > yU + /i~ with a branching ratio larger than the SM prediction does not 
necessarily mean a New Physics (NP) discovery. In such a case, the compatibility with the 



SM prediction has to be computed. Fig. 21 is the equivalent of Fig. 20 but with the SM 



rate for BR(£> S — > fx being considered as a background, and the signal corresponds to 
the NP contribution to BR(£> S — > We can see that for the same luminosity needed 

for a 3<t evidence of a SM-like signal, the LHC could alternatively claim NP at 3a if the NP 
contribution is of the order of 4 — 5 x 10~ 9 , i.e, if the actual BR(5 S — > p + ^) is O(8xl0~ 9 ). 
Finally, with the current uncertainties in fd/f s (7.9%) and in the SM prediction (8%), only 
values of BR(£> S — > that are at least 33% (55)% larger than the SM prediction can allow 

exclusion of a SM-like rate at 3 (5) a. 
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Figure 20: Required luminosity in order to provide a 3a evidence (orange) or a 5a discovery 
(green) of a given BR(5 S — > fi + n~) for LHCb and CMS combined. The luminosity is expressed 
in terms of the luminosity used in 15 , (0.34 fb -1 for LHCb and 1.14 fb _1 for CMS). 




Figure 21: Required luminosity in order to provide a 3a evidence (orange) or a 5a discovery 
(green) of a given NP contribution to BR(I? S — > /U + /i~) for LHCb and CMS combined. The 
luminosity is expressed in terms of the luminosity used for Ref. [15], (0.34 fb -1 for LHCb and 
1.14 fb" 1 for CMS). 



6 Conclusions 

The decay B s — > is known to be a very effective probe of SUSY models with large (> 30) 

tan/3, and its importance has been emphasised in numerous studies over the past decade. Due 
to its distinct signature, this decay can be searched for by three LHC collaborations: LHCb, 
CMS and ATLAS. Recently, searches by LHCb and CMS have been released, and have improved 
the upper limit on its branching ratio to BK(B S — > fi + n~) < 1.1 x 10~ 8 . Using this new bound, 
we performed a study of the constraints on the parameter space of five distinct SUSY models. 
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We emphasised that such indirect constraints can be stronger than those which are obtained 
from the ongoing direct searches for SUSY particles at the LHC For instance, in the CMSSM 
for tan /3 ~ 50, the SUSY particles have to be very heavy and in particular squarks cannot be 
lighter than ~ 1.2 — 2 TeV in order to be compatible with the upper limit on BK(B S — > fi + [i~). 
Nevertheless, in the scenarios we investigated here, in spite of the severe constraints we obtained, 
there is still room for SUSY contributions in large parts of the parameter space, especially for 
small tan (3. 

In addition, we considered an alternative observable which includes BR,(B S — > fi + fi~), 
namely a double ratio formed from the decays B s — > /i + pT , B u — > tu, D — >■ jiu and D s — > \iv j ru. 
The magnitude of the double ratio depends on the CKM matrix element \Vub\-, a parameter for 
which there is already considerable experimental information, and the prospects for further pre- 
cision in measurements of \V u b\ are promising. In contrast, the magnitude of BR(i? s — > /i + /i~) 
depends on the absolute value of the decay constant fs 3 , and thus a comparative study of the 
constraints obtained from these two observables is instructive. We showed that the double ratio 
can provide stronger constraints on the SUSY parameter space, and we advocate its use when 
discussing the impact of BR(£> S — > ji + //~) alone on SUSY models. 

The final integrated luminosity of the operation of the LHC at -y/s = 7 TeV is likely to 
be significantly larger than the amount that was anticipated at the start of the run. Both 
CMS and LHCb will have a chance to obtain a significant signal by the end of the run, even 
if BR(5 S — > is as small as the prediction in the SM. Throughout the run at y/s = 7 

TeV, the ongoing searches for B s — > will continue to compete with the direct searches 

for SUSY particles as a probe of the parameter space of SUSY models. 
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